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Some well-established examples of itinerant-electron ferromagnetism in one dimension occur in a 
Mott-insulating phase. We examine the consequences of doping a ferromagnetic insulator and cou¬ 
pling magnons to gapless charge fluctuations. Using a bosonization scheme for strongly interacting 
electrons, we derive an effective field theory for the magnon-holon interaction. When the magnon 
momentum matches the Fermi momentum of the holons, the backscattering of the magnon at low 
energies gives rise to a Kondo effect of a pseudospin defined from the chirality degree of freedom 
(right- or left-moving particles). The crossover between weak-coupling and strong-coupling fixed 
points of the effective mobile-impurity model is then investigated using a numerical renormalization 
group approach. 


PACS numbers: 71.10.Pm, 72.15.Qm, 75.10.Lp 

I. INTRODUCTION 

Ferromagnetism remains a challenging problem in 
physics, despite having been investigated ever since (and 
even before) the advent of quantum mechanics. It was 
Heisenberg who first realized that the phenomenon re¬ 
sults from an interplay between electron-electron interac¬ 
tions and the Pauli exclusion principle: 1 when two spin- 
polarized electrons occupy two orthogonal orbitals, their 
wave function must be spatially antisymmetric and van¬ 
ishes when they occupy the same position. This leads 
to a lower expectation value of the Coulomb repulsion 
for electrons with the same spin. The dependence of en¬ 
ergy levels on the relative spin orientation is often cast 
as an effective exchange interaction, which in this case is 
ferromagnetic as it favors parallel alignment of spins. 

The mechanism of direct exchange can be generalized 
to the many-electron problem. However, its implications 
for the ground state are not so clear. On the one hand, 
mean-field arguments predict that a gas of itinerant elec¬ 
trons with local repulsive interactions will spontaneously 
break spin rotational symmetry and become ferromag¬ 
netic for sufficiently strong interaction. The precise con¬ 
dition for the transition is given by the Stoner criterion: 
Up(Ep ) > 1, where U is the interaction strength and 
p(Ep) is the density of states at the Fermi level.® 

On the other hand, the Stoner criterion is not entirely 
reliable because the putative transition occurs in a non- 
perturbative regime of large U . Experimentally, the cri¬ 
terion remains controversial. In transition metal oxides, 
the existence of a ferromagnetic phase depends on de¬ 
tailed information about the band structure. 2 -! Recent ob¬ 
servations in cold-atom systems have concluded first in 
favor of and later against an interaction-driven ferromag¬ 
netic transition. 3 ^ 

Itinerant ferromagnetism has proved hard to estab¬ 
lish in microscopic models beyond mean-field approxima¬ 
tions. For many years, explicit proofs relied on peculiar 
conditions such as the limit of vanishing hole doping 
or the presence of flat bandsP In the domain of one¬ 


dimensional systems, there are even more constraints: a 
theorem due to Lieb and Mattis 7 rules out a ferromag¬ 
netic ground state for a number of models, including the 
paradigmatic Hubbard model. 

The Lieb-Mattis theorem does not hold for models with 
hopping beyond nearest neighbors. In fact, TasakP pro¬ 
posed a two-band one-dimensional model whose ground 
state can be shown to be fully polarized for a wide range 
of hopping parameters and finite repulsion. Remarkably, 
the proof relies on the condition of quarter-filling, for 
which the model was conjectured to be a Mott insulator. 
One might then wonder whether ferromagnetic order sur¬ 
vives in the metallic phase, reached by electron or hole 
doping. It is believed that it does, but the evidence relies 
mostly on variational methods and exact diagonalization 
for small chains. 2 9 There is stronger evidence based on 
density matrix renormalization group (DMRG) results 
for the single-band Hubbard model with next-nearest- 
neighbor hopping. 10 Only in the infinite-repulsion limit 
has metallic ferromagnetism been rigorously established, 
for a multi-band model with no flat bands. 11 

The nature of the ferromagnetic transition in one 
dimension has also been studied. Paramagnetic one¬ 
dimensional metals behave as Luttinger liquids, in which 
charge and spin excitations are described by two in¬ 
dependent charge and spin bosonic fields.®^ As the in¬ 
teraction increases and transition to a ferromagnetic 
phase sup posedl y occurs, the spin velocity must become 
negative. 10 13 14 Theories have been proposed to describe 
second-order transitions for Ising 13 15 and SU (2 f 1 ' 6 1 sym¬ 
metries. First-order transitions are also a possibility.® 7 ® 

In this work, we investigate the stability of doped 
Mott-insulating ferromagnets in one dimension at zero 
temperature by considering the creation of a magnon and 
its interaction with gapless charge fluctuations. This ap¬ 
proach follows the spirit of current experiments desi gned 
to investigate phase transitions in cold atomic gases!121^51 

We organize our presentation as follows. Our starting 
point in Sec. [H] is a generalized version of the models 
proposed in Refs. 3 and [9j In the weakly interacting 








regime, we bosonize the model and show that the charge 
sector is a Mott insulator for quarter-filling while the 
Luttinger liquid in the spin sector remains stable. Sec¬ 
tion |III| considers the model for the ferromagnetic phase 
in the strongly interacti ng reg ime. We use an alterna¬ 
tive bosonization scheme 21 22 J to identify the spin and 
charge excitations in the large -V limit as magnons and 
fermionic holons, respectively. In Sec. |IV| we analyze the 
low-doping limit of the metallic phase using an effective 
field theory with magnon-holon interactions dictated by 
symmetry considerations. A similar approach has been 
used to study magnons in a spinor Bose liquid. 23 We 
find that, when the magnon momentum is commensurate 
with the holon Fermi surface, the scattering between op¬ 
posite Fermi points gives rise to infrared singularities akin 
to the Kondo effect “This phenomenon has been previ¬ 
ously discussed in the context of a mobile impurity in a 
Luttinger liquid.^ 26 In Sec. [v| we proceed to studying 
the low-energy fixed points of our eff ective model using 
the numerical renormalization group. 27 28 Finally, Sec¬ 
tion ED summarizes our results. 


II. WEAKLY INTERACTING REGIME 


In Sec. [I] we briefly mentioned Tasaki’s model, a one¬ 
dimensional model for which ferromagnetism has been 
rigorously established.^ Tasaki’s Hamiltonian is 

H = E ^' c J,<7 c j> + u E t n u> (!) 

i,j,cr i 


where is the annihilation operator for an electron 

with spin a =t,4- at site h n i,o- = C I a c i a * s the number 
operator, U > 0 is the on-site repulsion, and the hopping 
amplitudes Uj are defined as follows: 


U,i+1 = ^ 5 

U,i+ 2 = ^ 

U,i+ 2 S, 

ti,i = v, 

ti,j — P 


Vi, 

if i is even, 
if i is odd, 
if i is odd, 
otherwise. 



Figure 1. (color online) Tasaki’s one-dimensional model. 
The hopping between nearest neighbors is t'. The hopping 
between next-nearest neighbors is t (— s) for the even (odd) 
sublattice. The odd sublattice has chemical potential V. The 
on-site electronic repulsion in both sublattices is U. 


Let us describe the low-energy excitations of Tasaki’s 
model, starting with the paramagnetic phase at weak 
coupling. As a first step, we diagonalize the Hamilto¬ 
nian for U = 0, which we denote by Hq. Using the mo¬ 
mentum representation defined on the even and the odd 
sublattices 

%-^E^W, ( 2 ) 

C2j + l,a = y= E etk(23 + 1)b k,a, (3) 


with k G [—7t/2,7t/ 2] in the first Brillouin zone, we can 
rewrite Hq as 


H 0 = E ^tcos(2k)al a ci ktCr - 2scos{2k)b\ a b k ^ 

k,cr 

+2 1' cos k{b\ a a k ^ + H. c.) . 


(4) 


We see now that Tasaki’s model is a particular case of a 
model of the form 


k,cr 

+a(k)(bl a a k)<T + h.c.) . 


(5) 


This model can be diagonalized through a rotation 


a /c,cr — COS 

b k ,a = sin 



9+,k,a ~ sin 


g+,k ,<j + cos 



(6) 

(7) 


Fig- [T] illustrates the Tasaki lattice. For i / -s, the 
lattice has two sites per unit cell. Tasaki demonstrated 
the existence of a ferromagnetic phase at quarter-filling 
when t and U are sufficiently large compared to s and for 
the particular choices V = 0 and t' = y/2 (t + s). Under 
these conditions, it can be shown that the expectation 
value of the Hamiltonian in an arbitrary state satisfies 
the inequality (H) > —2 (s + t) N , where 2 N is the total 
number of sites. The lower bound happens to be pre¬ 
cisely the energy of the fully polarized state, which proves 
that the latter is one of the ground states P The ground 
state is (2 S to t + l)-fold degenerate, with S to t = AT/2 at 
quarter-filling, due to spontaneous symmetry breaking of 
the SU( 2) symmetry. 


where 


2 a(k) 

tan Ok = 


e«(fc) - e&(fc)' 


(8) 


We obtain two bands 

Fo = E [ e +( k )9+,k,*9+,k,<r + e-{k)gl ik:(7 g-, k ,a , (9) 


k,a 

with dispersion relations 

£o(fc) + £b(ft) , 


e±(k) 


i / /0/2 / b>\ 

~e a (k) - e b (k)~ 


2 


-I 2 


( 10 ) 














The noninteracting ground state is paramagnetic. At 
quarter-filling the single-particle states in the lower band 
with \k\ < kp = 7r/4 are occupied (corresponding to 
half-filling of the lower band). The neutral excitations 
are electron-hole pairs on top of this ground state. 

Now we consider the weakly interacting regime 


U < vf, where vp 


k=k f 


is the Fermi ve¬ 


locity. The weak-coupling condition justifies lineariz¬ 
ing the low -ener gy spectrum around ±kp. Through 
bosonization, 29 30 one can map the linearized version of 
the kinetic energy Hq onto an effective Hamiltonian 


Kl 


E 


/ 


dx Y i 


n 2 + (d x </> v y 


(ii) 


where (j) c and n c (</> s and U s ) are the charge (spin) 
canonically-conjugated bosonic fields. It can be shown 
that d x (j) c {x) ~ p(x) and d x (f) s (x) ~ ra(x), where p(x) is 
the charge density and m(x) is the local magnetization. 
The remarkable feature of bosonization is that, once the 
interacting term is bosonized as well, charge and spin 
fields remain noninteracting, the Hamiltonian of Eq. (11) 
being only slightly modified® 


H ll = 


7 y- r> o ** 


dx Y 


K„nl + T (, d x ^) 2 


( 12 ) 


Hamiltonian Hpp is known as the Luttinger model; v c 
and v s are the charge and spin velocities; and K c and 
K s are the Luttinger parameters in the charge and spin 
sectors, respectively. In the noninteracting case, v c = 
v s = vp and K c = K s = 1. For repulsive interactions, we 
have K c < 1, while spin SU(2) symmetry fixes K s = l. 31 
The interaction has the effect of modifying the velocities 
and Luttinger parameters, leading to the phenomenon 
of spin-charge separation: in the low-energy limit, the 
elementary excitations are independent charge and spin 
collective modes. 

At weak coupling, we can determine the U dependence 
of charge and spin velocities by directly bosonizing the 
interaction term in Eq. 0- For this purpose, we first 
promote the site operators to fields. Defining a field is not 
straightforward as the expressions for c^j a and i, a 
are different. Combining Eqs. 0, 0, 0 and Q to 
express the site operators in terms of g±,k,a yields 


C 2 j,< 


v-^ e 


ik(2j) 


0 * 


~-2.y7^ sli yy ) 9 -,k,a + -, (is) 


e ik(2j+l) 

x -—— cos — 

y/N V 2 


C 2 j+l,a ~ 2^ - 7^7 - COS ( — ) 9-’k,° + ’ • • J ( 14 ) 


where we have omitted terms proportional to g+,kp be¬ 
cause the latter involve high-energy excitations, in the 
upper band. 


We can combine Eqs. (13) and (14) into a single field 
as follows 


ip<r(x) ~^2e lkx f(k,x)g- tkt0 


(15) 


where we have introduced the function 

l-he i7rx fe k \ l- e i7rx 
cos 


f(k,x) = 


sm I % ) , (16) 


which reduces to f(k,x) = cos(0^/2) [f(k,x) = 

sin(0fe/2)] for the positions x belonging to the even (odd) 
sublattice. 

The next step is to expand the electron field into right 
and left movers 


4>a(x) ~ f(k F ,x) [e lkFX Tp L ,a(x) + e lkFX 'ipR !<T (x)] , 


where 


j(x) 


1 


+oo 


„±ipx 

e 9 (R,L),p,ct- 


(17) 


(18) 


Here g(R t L),p,a are the annihilation operators of the lin¬ 
earized branches around the right and left Fermi points. 

In the continuum limit, the interaction term in Eq. 0 
becomes UJ2j n j,t n j,i ~ U J dxip^ip^ip^ip^. The proce- 
dure to bosonize the interaction is now almost identical 
with that for the Hubbard modelp^ except for the alter¬ 
nation between even and odd sublattices that introduces 
the factor of f(kp,x). From now on, we follow the same 
steps as for the Hubbard model, but must be careful to 
take into account the oscillations of f(kp,x). 

The uniform part of the density operator, combined 
with the x-independent part of f(kp,x ), gives rise to 
terms in the effective Hamiltonian which are quadratic 
in the bosonic fields, as in Eq. ( [I2| ). Compared with the 
result for the Hubbard model, the effective U is renor¬ 
malized into an effective U = \ [5 + cos( 26 > / Cf )] U. (The 
proportionality factor actually depends on the choice of 
/(fc, x) introduced in Eq. (16), but the important result is 
that U ~U.) Therefore, one must just replace U by U in 
the expressions for the Luttinger parameters of the Hub¬ 
bard model. In particular, the spin velocity for U vp 
is known to be 31 


i--?0 


1/2 ( u \ 


1/2 


v s « Vp [ l-l Vp 1-. (19) 

TTVp J \ TTVp I 

The Luttinger liquid phase becomes unstable when 
v s —» 0, which we can interpret as a sign of a phase 
transition to a state with spontaneous magnetization. 13 
In principle, one could search for this instability by ex¬ 
trapolating the result in Eq. (19). However, the condition 
U = TTVp is not compatible with U <C vp, required for 
weak-coupling bosonization to hold. In the next section 
we shall consider the effective field theory in the strong- 
coupling limit. Before doing so, we now argue that the 
model is a Mott insulator at quarter-filling, even at arbi¬ 
trarily small U. 

The insulating behavior is due to Umklapp scat¬ 
tering, which becomes commensurate in the two- 
sublattice system when kp = 7r/4. The Umk¬ 

lapp operator stems from the interaction term ^ 


















U f dx[f(kp, x)] 4 e z4kFX ^ R ^ R ^ R c.. The os¬ 
cillating component of f{kp,x) can be combined with 
the oscillations of the fermionic field to generate terms 
proportional to Often one argues that these 

terms oscillate rapidly and can be neglected in the low- 
energy Hamiltonian. However, precisely for quarter- 
filling, the oscillations cease and these terms must be 
kept. The bosonized version of the nonoscillating Umk- 
lapp process is 


HlJmklapp — 9 1 dx cos ^\/ 87 ri 


( 20 ) 


where g is the coupling constant. From the 
renor maliza tion-group analysis of the sine-Gordon 
modelC^ES] ^ is known that this operator has scaling 
dimension 2 K c and is relevant for arbitrarily weak re¬ 
pulsive interactions. Its effect is to open up a gap 
A c in the charge sector, which for small q scales as 

A particularly simple result is obtained at the Luther- 
Emery point K c = 1/2, in which the sine-Gordon 
model can be refermionized into noninteracting spinless 
fermions . 34 The gap can be seen explicitly in the mas¬ 
sive relativistic dispersion 


<T,± 


(p) = ivVcP ) 2 + A 2 . 


( 21 ) 


The positive (negative) energy, as measured from the 
chemical potential, refers to a completely empty (filled) 
band of free fermions which are the elementary excita¬ 
tions in the charge sector. The gap between the bands 
at the Luther-Emery point is 2A C ^ \g\. 

In our case, a caveat is necessary. Setting kp = 7 r /4 
in Eq. (16) yields g = 0 since 0k =7r /4 = ^sgn(£') for 
Tasaki’s model (see Eq. <©)• However, this bosonization 
procedure only predicts the coupling constants correctly 
to first order in U. Since the Umklapp operator is not for¬ 
bidden by any symmetries in Tasaki’s model, we expect 
it to be generated at higher order, most likely g ~ 0{U 2 ). 
Therefore, the effective field theory predicts the system 
to be paramagnetic at weak coupling, i. e., a Luttinger 
liquid in the spin sector, and to become a charge insulator 
at quarter filling. 


III. STRONGLY INTERACTING REGIME 

We have seen in Sec. [H] that there is no ferromag¬ 
netic transition for the model of Eq. © at small U . On 
the other hand, at least for Tasaki’s model at quarter¬ 
filling we know that the ground state becomes a fully 
polarized ferromagnet at sufficiently large U P Since we 
expect the charge gap to increase monotonically with 
U , there should be a transition from a paramagnetic 
Mott insulator (with gapless spin excitations described 
by the Luttinger model) to a ferromagnetic Mott insula¬ 
tor (with gapless magnons due to the broken SU( 2 ) sym¬ 
metry). Without discussing the nature of the transition 


(whether first or second ordeiE^HMII]), we now consider 
the strongly interacting regime and take the existence of 
a ferromagnetic ground state for granted. We shall start 
from the insulating phase at quarter-filling, but also in¬ 
vestigate the consequences of doping into the metallic 
phase. 

The difficulty in treating the problem in the strongly 
interacting regime is that standard bosonization is 
not applicable. Nonetheless, an alternative bosoniza¬ 
tion scheme for strongly interacting electrons has been 
developed . 21 This approach starts from the picture that, 
in the limit of infinite repulsion, electrons cannot move 
past each other; as a result, the spin of each electron 
is confined and therefore frozen. One then writes down 
an effective model for spinless fermions ( holons ) in the 
charge sector. Including corrections due to forward scat¬ 
tering at large, finite U gives rise to an exchange inter¬ 
action J, which allows one to treat the spin degrees of 
freedom as an effective spin chain. In this scenario, spin- 
charge separation still holds . 21 

We follow this strategy, but adapt it to the problem 
considered here. The main difference is that, in the pres¬ 
ence of next-nearest-neighbor hopping, an electron can 
hop around another electron that occupies a nearest- 
neighbor site. Therefore, the spin is not confined even 
for infinite repulsion. To apply the picture of Ref. 21 ] we 
consider the more general Hamiltonian in Eq. § with 
chemical potential V ^ 0. Recall that Tasaki’s proof 
of a ferromagnetic ground state only applies to the case 
V as 0, t' as y/2(t + s). However, adding the staggered 
chemical potential does not change the symmetry of the 
Hamiltonian, and it has been argued that the ferromag¬ 
netic phase is observed also for V ^ 0. 2 9 

Here, we take V to be negative and large, \V\ £, s, 

to constrain the electrons to move within the odd sub¬ 
lattice, with hopping amplitude —s. Together with a 
strong on-site repulsion, U s, this condition suppresses 
exchange processes, freezing the electron spin degree of 
freedom in the limit |H|,L/ oc. The charge sector 
can then be described by spinless fermions, which we call 
holons from now on. In this regime it is easy to see that 
the system is a Mott insulator at quarter filling, since it 
corresponds to half-filling of the odd sublattice, with a 
gap in the charge excitation spectrum of order U s. 
At the same time, we can think that these gapped holons 
descend from the fermions for the sine-Gordon model at 
the Luther-Emery point discussed in Sec. [TTJ Extending 
the dispersion in Eq. (21) to the large-gap regime and 
expanding for \p\ A c ~ U, we have 


e c,±(p) ~ ±(A C + 7 ±p 2 ). 


( 22 ) 


Here and are parameters inversely proportional to 
the holon mass in the upper and lower Hubbard bands, 
respectively. In contrast with Eq. (21), we allow for 


7 + 7 ^ 7 _ because the effective model in the large -U limit 
has no Lorentz invariance nor particle-hole symmetry. In 
fact, we expect 7 _ ^ s and 77 ~ £ from hopping of holes 
in the odd sublattice and particles in the even sublattice. 












Figure 2. (color online) Spin exchange mechanisms for the 
Tasaki lattice. The first mechanism is the same that takes 
place in the Hubbard model and is always antiferromagnetic. 
The spin flip occurs through an intermediate state of addi¬ 
tional energy U 1 as represented in the upper box. In the 
second mechanism, represented in the lower box, an electron 
moves to the other sub-lattice, creating a vacancy in the first 
sub-lattice for another electron. Next, the second electron 
takes that vacancy. Finally, the first electron occupy the site 
of the second one. The exchange constant is ferromagnetic if 
—s < 0 . 


Now we take into account that the repulsion and 
the chemical potential are finite and allow spin fluctu¬ 
ations. Treating virtual hopping processes by perturba¬ 
tion theory leads to an effective spin exchange interac¬ 
tion, as illustrated in Fig. [2] The resulting Hamiltonian 
is H s = JJN • Si+i, where is the spin operator of 
the i-th electron in the chain, and the exchange constant 
is given b} 09] 


J = 4 



st^\ 

V 2 ) ' 


(23) 


Note that J can be ferromagnetic only if s > 0 (recall 
that the hopping amplitude in the odd sublattice is — 5 , 
hence negative). We will take J < 0 from now on. In 
this case, the ground state of the effective spin chain is a 
fully polarized ferromagnet. 

Taking the ground state to be fully polarized along 
the £ direction, the spin Hamiltonian can be mapped 
into magnon excitations through the Holstein-Primakoff 
transformation!^ 


S^S-Bfa, (24) 

SJ =B]{2S-B]B,)' / ^V2SB], (25) 

S+=( 2 S-B]B j )‘ / 2 B j! »^B j , (26) 


where S is the spin quantum number (S = 1/2 for elec¬ 
trons), and the magnon operators Bj obey a bosonic al¬ 
gebra. We see from Eq. (24) that the number of magnons 
is directly related to the magnetization. 


Within standard linear spin-wave theory, the spin 
Hamiltonian describes noninteracting magnons and takes 
the form 


(27) 

k 

with uo(k) = 2 JS (1 — cos k) ~ JSk 2 for k < Cl. The 
dispersion is gapless because the magnon is the Gold- 
stone mode of the spontaneously broken SU( 2 ) sym¬ 
metry. From now on, we write simply uj(k) ~ A k 2 at 
low energies. More generally, the effective A away from 
the strong-coupling limit can be determined within a 
random phase approximation for the spin-spin correla¬ 
tion function . 37 Neglecting interactions between magnons 
and gapped charge fluctuations, the ferromagnetic Mott- 
insulating phase at quarter-filling is stable as long as 
A > 0 . 

We are interested in the effects of coupling magnons 
to gapless charge excitations which arise when we move 
away from quarter-filling. Since the density of holons is 
directly related to the density of electrons, adding elec¬ 
trons imply adding holons to the upper Hubbard band. 
The configurations before and after doping are illustrated 
in Fig. [3] The doping introduces two Fermi points at 
=b pf, where pp is related to the average charge density 
P = E j,A n i,a) by Pf = tt(/ 9- 1/2). At energy scales be- 
low the charge gap, we can ignore the lower band. More¬ 
over, in the low-doping regime (described by the the¬ 
ory of the commensurate-incommensurate transition 35 ) 
we can neglect holon-holon interactions, which are weak 
even when K c 7 ^ 1/2. The charge Hamiltonian is then 
simply 


H r 


^2 e c(p) c t c P’ 

p 


(28) 


where cj, creates holons in the upper band and e c (p) « 
7 p 2 (omitting the constant A c in the energy), with 7 = 
7 _i_. (Likewise, removing electrons would create holes in 
the lower band. In this case one would have to replace 7 + 
by 7 _ for the low-energy excitations.) We recall that in 
the strongly interacting regime we expect 7 ^ £, whereas 
A ~ J with J given in Eq. (23). Therefore, our analysis 


is valid in the regime A 7 , i.e. when the magnon mass 
(in the Galilean sense of a quadratic dispersion) is larger 
than the holon mass. 

We now discuss the interaction between magnons and 
holons. The form of the interaction can be guessed from 
symmetry, but before doing so we propose a physical 
mechanism from which it can be derived. In general, 
the exchange constant J is a function of the average 
charge density, J = J(p). Following Refs. 38l and l39l 
we consider that in the metallic phase the exchange in¬ 
teraction between neighboring spins is a function of the 
local density: H s = YU J(^)S* • S* + i. In our case, 
pi can fluctuate due to the motion of the dilute elec¬ 
tron gas in the even sublattice. Next, we expand J(pi) 
around the average density before doping (po = 1 / 2 ): 








Before doping 


After doping 


p-q 




Figure 3. (color online) Holon bands before and after doping 
the ferromagnetic insulator. At low energies, the dispersion 
is approximately parabolic. The doping introduces two Fermi 
points, — pf and pf- 


P 



k + q 
k 


C (k,q) ~ k(k + q) 


Figure 4. Magnon-holon interaction. Solid (dashed) lines 
represent holons (magnons). The momentum dependence of 
the coupling function £ is dictated by symmetry. 




Figure 5. First- and second-order diagrams in the magnon 
self-energy due to magnon-holon interaction. 


J(pi) ~ J(po) + jp~\ p ._ po $Pi- The term proportional to 
the derivate of J leads to the spin-charge interaction 


Ho = Y, 


SJ(pi 


Spi 


8pi&i • Si 


;+i- 


(29) 


pi=po 


Replacing 8pi by the holon density operator of the upper 
band, we can write 


h cs — 4 c i s * ■ s 


t+i> 


(30) 


where a = — ^ 


p o 


and Ci is the Fourier transform of c v in 


Eq. (28). Using the Holstein-Primakoff transformation 


again and switching to the momentum representation, 
we obtain the quartic magnon-holon interaction Hi = 
T,k, q , P C{k,p,q)cl_ q c P Bl +q B k , with 


C (k,p,q) = -a 


e ik _|_ e ~i{h+q) 


ak(k + q). 


(31) 


Note that the coupling function is proportional to the mo¬ 
menta of the incoming and the outgoing magnons. Next 
we discuss how the same form is enforced by symmetry. 

Due to spin SU( 2) symmetry, the total magnetiza¬ 
tion S* ot = J2j Sj * s a conserved quantity in our model. 


Therefore, according to Eq. (24), the number of magnons 


must also be conserved. This implies that the magnon- 
holon interaction vertex £ must contain one incoming and 
one outgoing magnon legs (in contrast with the more of¬ 
ten encountered coupling between electrons and phonons 
or photons). By charge conservation, there is also one 
incoming and one outgoing holon. 

Let p and k be the momenta of the incoming holon and 
magnon, respectively, and p — q and k + q the momenta 
of the outgoing holon and magnon. In principle, the cou¬ 
pling function ( can depend on p, k and q. Parity symme¬ 
try requires that the expansion of ( for small p,k,q con¬ 
tains only even powers of momenta. Furthermore, there 


should be no coupling to magnons with zero momentum 
(Goldstone bosons) since their presence only amounts to 
a uniform rotation of the magnetization. Thus, £ —> 0 if 
k-^Oork + q^O. This rules out a constant term in 
The lowest-order term is then a product of two momenta, 
so we must have £ oc k(k+q). Therefore, the generic form 
of the magnon-holon interaction in the long-wavelength 
limit is indeed 


Hi = a J2 k (k + q)c ] p _ q CpBl +q B k . (32) 

k,p,q 


This result is summarized in Fig. [4j 

Our complete Hamiltonian is H = H c + H s + Hj. To 
find out how the interaction affects the magnon spec¬ 
trum, we calculate the magnon self-energy at first and 
second order in a. The corresponding diagrams are illus¬ 
trated in Fig. [ 5 ] In the first-order (“tadpole”) diagram, 
the holon loop is proportional to the holon density and 
the interaction vertex contributes with a factor of k 2 . 
Therefore, the first order diagram only renormalizes the 
magnon mass, with 8 A ~ 0(a). 

The second-order diagram is more interesting. To 
make the integrals well defined, we introduce a momen¬ 
tum cutoff D around the Fermi points — pp and pp for 
magnons and holons. We divide the self-energy into two 
contributions, ^ 2 \k,uj) = (fc, cu) + (k, u). Both 

contributions contain logarithmic singularities in their 
real parts: 


5?4 2) - a 2 


(u) — A k 2 ) log ^ 


D 


uj — A k 2 
2 


~ a 2 [u - \(k - 2 p F ) z ] log 
+(Pf -Pf)- 


D 


uj — A (k — 2 pf) 


(33) 


The first contribution arises from low-momentum scatter¬ 
ing and is typical of the orthogonality catastrophe, which 


























has already been studied in the case of bosonsThe 
second contribution corresponds to a momentum trans¬ 
fer of 2 p F \ it diverges on-shell only if \k db 2 p F \ ~ |fc|, 
where k is the magnon momentum. These singularities 
tell us some important mechanism is taking place when 
magnons and holons scatter around the Fermi points. 


IV. EFFECTIVE FIELD THEORY FOR 
MAGNON-HOLON INTERACTION 

In this section, we derive an effective field theory for 
the scattering between holons and magnons with momen¬ 
tum close to the Fermi points ±p F . Our goal is to en¬ 
lighten the mechanism behind the infrared singularities 
encountered in Sec. m and to set the stage for going 
beyond perturbation theory using the numerical renor¬ 
malization group in Sec. [V| 


A. Chirality Kondo effect 


We start by rewriting the free Hamiltonians of Eq. (27) 


and ( [28] ) in the continuum limit in terms of holon and 
magnon fields: 


H c = 7 J dx d x ip\x)d x ip(x), (34) 

H s = X J dx d x B\x)d x B{x). (35) 

where ip(x) = Fp c P eipx and B ( x ) = 7 ^ J 2 k B ke ikx , 
with L the system size. 

We now restrict to low energies compared with the 
holon Fermi energy \p 2 F . In this regime holons can only 
be scattered in the vicinity of the Fermi points. We ex¬ 
pand the holon fields in terms of left movers and right 
movers, 


« e~ ipFX ^ L {x) + e ipFX ^ R {x). (36) 

We also focus on magnon states with momentum near 
±Pi?, and expand 

B(x) « e~ ipFX B L (x) + e ipFX B R (x). (37) 

We define two-component spinor fields which combine 
right and left movers 

*«-($:!)• nw =(£«)• (38) 



Figure 6. (color online) Different kinds of fermion-magnon 
scattering around the Fermi levels. The indices L and R indi¬ 
cate whether the particle momentum is closest to — p F or p F . 
By defining a chirality pseudospin such that R = 7 and L = 4, 
the scattering can be identified with pseudospin exchange in¬ 
teraction (as represented by the diagrams on the bottom). 
The solid line represent holons and the dashed line represent 
magnons. There are in fact six diagrams, but the remain¬ 
ing ones can be obtained from the ones presented through a 
parity transformation. 


where v = 2 ^p F and u = 2A p F are the holon and magnon 
velocities, respectively, and r z is the z Pauli matrix in 
the internal chirality space. As discussed in Sec. |III[ we 
expect A 7 in the limit J s, t. Thus, we are in a 
“slow magnon” regime u <C v. The spinor representation 
suggests regarding the chirality indices as pseudospins, 
R = t and L = l (not to be confused with the original 
electron spin a). We will see soon that this picture is 
useful to interpret the interaction, but it is important to 
keep in mind that the physical meaning of R (L) is that 
the particle carries momentum close to p F (—p F ). 

Now we must propose the form of the magnon-holon 
interaction. As before, we appeal to symmetries. As 
we discussed in Sec. |III[ conservation of magnetization 
implies that a local interacting term must have the gen¬ 
eral form 'ipl i (x)'ip r 2 (x)Bl 3 (x)B r 4 (x) in order to conserve 
the number of magnons. In this notation, r represents a 
pseudospin index (R or L). Momentum conservation im¬ 
pose constraints on the r’s: only combinations such that 
ri — 7*2 + rs — rq = 0 are allowed. These are: LLLL and 
RRRR\ LLRR and RRLL\ and LRRL and RLLR. We 
split them into pairs related by a parity transformation 
(where L R). To ensure parity symmetry, coupling 
constants for each element of a pair must be the same. 
The corresponding scattering processes are illustrated in 
Fig. P 

With these restrictions, it is not difficult to show that, 
in the spinor representation, the interacting part of the 
Hamiltonian assumes the form 


= y J dx 7 / dx 

+ J dx (4/^r + 4/ + H. c.), (41) 


Linearizing both holon and magnon dispersion for k,p ~ 
=b p F (and measuring the holon dispersion from the Fermi 
energy), we can write 


H c 

H s 


J dx T^(x) (—ivr z d x ) T(x), 

J dx Cfl(x) [oj(pf) ~ iur z d x \ Q(x), 


(39) 

(40) 


where V, J \\, J_\_ are three independent coupling constants. 
Since this interaction term must arise from linearizing 
Eq. (52), we expect the coupling constants to be of the 


order of 


(XPp. 


In terms of pseudospins, the first term in Eq. (41) 
describes the potential scattering between holons and 
magnons, while the remaining terms describe an 
anisotropic exchange interaction between holon and 







magnon pseudospins. Recognizing S p(x) = ^(x)^(x) 
and S m{x) = (x)^Q(x) as pseudospin vector densities, 

we further rewrite Hj as 


Hr = V J dx + J|| J dx S z F (x)S z M (x) 

+ ' J j J dx [Sf(*)Sm(x) + H. c.] , (42) 


where the label F is for (fermionic) holons and the label 
M is for magnons. 

The Hamiltonian in Eq. ( [42] ) can be viewed as a Kondo 
interaction^ between a finite density of fermions (repre¬ 
sented by the holons) and a number of mobile impurities 
(represented by the magnons). The connection with the 
Kondo model is made more clear in the limit of vanishing 
magnon density (which is the interesting regime to in¬ 
vestigate the stability of the ferromagnetic state against 
infinitesimal perturbations). Since magnetization is con¬ 
served, we may restrict our problem to the subspace with 
a single magnon. In this case we can rewrite the Hamil¬ 
tonian using first quantization for magnon operators: 


H s = uj(p F ) + 2 uS z M P, 

H! = (*)*(*) + J l{ S z F (X)S z M 

, J± 


s$(x)s: 


M 


H. c. 


(43) 


(44) 


where X and P are the position and momentum oper¬ 
ators of the magnon, respectively, obeying [X,P] = i. 
Note that the Hamiltonian is translationally invariant, 
but the interaction is local at the position of the magnon. 
This kind of dynamics has been encountered in the cou¬ 
pling of a mobile impurity with a Luttinger liquid. 25 

Next, we apply a Galile an transformation to the 
magnon reference frame! 41 42 The transformation is 

JJ = e -iX f dx &(- id x ) V , (45) 


and is such that 

T(x) = U^(x)U = V(x - X), (46) 

P = U ] PU = P - J dx&i-idx)®. (47) 

As a result, the transformed interaction Hj = WHjU 
becomes 


Hj = (o)^(o) + j\\Sfr(o)S[ 

+ J Y [^(0)5m + H. c.] 


M 


(48) 


The magnon-holon interaction is now restricted to x = 0, 
since we are in a frame whose origin moves along with 
the magnon. For consistency, we must also transform the 
kinetic energy terms H c and H s . The free holon Hamil¬ 
tonian remains unchanged, i. e., H c = H c . However, the 
magnon term in Eq. (43) gets modified: 


H s = U ] H 8 U 

= uj(pf) + 2uS z M P + 2iuS z M J dx ^dx'b. (49) 


At this point we have eliminated the magnon position 
X from the transformed Hamiltonian H = H c + H s + Hi . 
Therefore, [IT, P] = 0, implying that H and P can be 
simultaneously diagonalized. This is a consequen ce of 
translational invariance of the original Hamiltonian. 4 dE21 
We shall focus on the subspace with eigenvalue P — 0, 
which corresponds to a magnon with energy exactly equal 
to In this case the first term in Eq. (49) vanishes. 

Deviations from the P = 0 condition amount to an extra 
term proportional to S Z M , which is equivalent to an ef¬ 
fective magnetic field acting on the magnon pseudospin. 
While we have got rid of X, the price we paid is the in¬ 
troduction of the second term in Eq. (49). This term 


couples the magnon pseudospin to the total momentum 
operator for holons. 

We can now add H s for P = 0 to the free holon Hamil¬ 
tonian H r : 


H r + H s — id 


( Pf) + J dxiphiv - 2uS z M )(-id x )tp f 


/ 


dx ipl(v + 2 uS z M )(id x )i> L . 


(50) 


We further rewrite the right-hand side of Eq. ([50]) by 
taking ^l(^) (this is equivalent to relabeling 

holon states with momentum p as — p for the left branch): 


H c -\-H s — id 


(Pf) + J dx^(v-2uS z M T z )(-id x )\l>. (51) 


We emphasize that r z in Eq. (51) acts in the holon pseu¬ 


dospin space. Recall that v is the holon velocity around 
the Fermi points ±p_p. Eq. (51) shows that the holon 


velocity in the transformed Hamiltonian is now given by 
v d= u depending on whether the holon and magnon have 
the same or opposite pseudospins, i. e. move in the same 
or opposite directions. 

Finally, dropping the constant ud(pp) and adding the 
transformed interaction of Eq. (48), we arrive at our effec¬ 


tive Kondo-type Hamiltonian for the scattering of a sin¬ 
gle magnon with momentum close to ±Pf by low-energy 
holons: 


+oo 

H= Y, uS z M T z )k<H k 

k— — oo 

+ ^[^(0)% + H. c.] 

+J\\S z f (0)S z m + V9'(PM0). (52) 


B. Perturbative renormalization group analysis 


If the magnon had zero velocity, i. e., u = 0, the model 
of Eq. (52) would correspond precisely to the Kondo 
model, which describes a Fermi sea of electrons cou¬ 
pled to a spin-1/2 impurity at the origin through an 
anisotropic exchange interact ion Here the model also 
includes the potential scattering V, which is generally 













present in the absence of particle-hole symmetry. It is 
well known that perturbation theory in the Kondo inter¬ 
action gives rise to infrared divergences like the ones we 
encountered in Eq. ( [33] ) . The Kondo model has been in¬ 
vestigated non-perturba tively using the numerical renor¬ 
malization group (NRG)P? 28 and solved exactly through 
Bethe ansatz.® From the renormalization group (RG) 
point of view, the logarithmic singularities that arise in 
perturbation theory can be recast in the renormalization 
of effective coupling constants at the appropriate energy 
scale. For u = 0, the exchange coupling constants flow to 
strong coupling if J± > — Jy. The low-energy fixed point 
of the model can be understood in terms of the forma¬ 
tion of a singlet between the localized impurity and one 
electron from the conduction band, which decouples from 
the remaining electrons. In other words, the fermion spin 
screens the impurity spin. 

In our case, the renormalization of the effective ex¬ 
change interactions decides the fate of the magnon in the 
low-energy limit. We first approach the problem by ap¬ 
plying a perturbative RG approach to interaction (41) 


by Eq. (56). By analogy with the original Kondo ef- 


with u > 0. The procedure is similar to Anderson’s poor 
man’s scaling,^ generalized to include a momentum cut¬ 
off in the magnon sub-band and to respect total momen¬ 
tum conservation in the virtual magnon-holon scattering 
process. We obtain the set of RG equations: 


_ vJ±(J\\ - V ) 

d£ 27 r(v 2 — u 2 ) 

d(J\\ - V) 


/ 2 


d£ 


d(J\\ + V) 


27 r(v + u) ’ 

J? 


d£ 2tt(v — u ) ’ 

du uJj_ 

d£ 87 t(v 2 — u 2 )' 


(53) 

(54) 

(55) 

(56) 


where d£ = dA/A denotes the infinitesimal reduction of 
the ultraviolet cutoff in the RG step. 

To analyze the RG flow, we first note that to second or¬ 
der in perturbation theory the renormalization of Jy + V 
is controlled by J± and does not feed back into the beta 
functions for J± and Jy — V. We can then focus on Eqs. 
(53) and ([54]) . In the slow magnon regime u < v, these 


coupled equations define a Kosterlitz-Thouless type flow 
diagram in the (Jj_, Jy — V) plane, with a flow to strong 
coupling for J± > — (Jy — V). Furthermore, Eq. (56) 
implies that the magnon velocity decreases with the RG 
flow. Here we should remark that we can follow the flow 
to arbitrarily low energies only in the case P = 0, i. e. 
when the magnon momentum exactly matches the holon 

IVA 


Fermi momentum pp. As discussed in Sec. 


devi¬ 


ations from this condition are equivalent to an effective 
magnetic field coupled to the magnon pseudospin, which 
cuts off the RG flow at the scale set by P. 

The enhancement of the effective couplings suggests 
that in the low-energy limit the magnon is strongly 
backscattered between states with momentum ±p_p, and 
its group velocity may eventually vanish as suggested 


feet, we expect that at the low-energy fixed point the 
magnon will form a “chirality singlet” with one holon, 
as represented by the state |$) = A= f dx[B R (x)'if R (x) — 

B\ j (x)'ip R (x)]\0). This singlet involves a pair of (distin¬ 
guishable) particles which are delocalized in space and 
move with opposite momenta — somewhat akin to a 
Cooper pair in BCS theory. Note that the particle-hole 
symmetry of our model is broken both by V 7 ^ 0 and 
u 7 ^ 0 ; thus, the putative low-energy fixed point should 
contain a marginal potential scattering operator that ac¬ 
counts for a non-universal phase shift in the remaining 
holon states. In Sec. [V]we shall investigate the nature of 
the low-energy fixed point using the NRG method. 

V. NUMERICAL RENORMALIZATION GROUP 


The numerical diagonalization of the Hamiltonian (52) 


follows the standard NRG method. The procedure hav¬ 
ing been thoroughly detailed , 27 28 only cursory descrip¬ 
tion is necessary, meant to prepare the discussion of the 
flow in renormalization-group space that constitutes the 
central object of this section. 


A. Procedure 


To diagonalize the Hamiltonian (52), we will rely on 


the numerical renormalization-group procedure, an ap¬ 
proximate iterative method that depends on strictly con¬ 
trollable approximations to yield a set of eigenvalues and 
eigenvectors from which essentially exact physical prop¬ 
erties. More important from our perspective, the it¬ 
erative diagonalization can be regarded as a sequence 
of renormalization-group transformations and therefore 
accurately describes the flow of the Hamiltonian in 
renormalization-group space. 

To define the renormalization-group transformation, 
we need a scaled truncated version of the Hamiltonian. 


Our initial goal is to project Eq. (52) on a finite basis. 


Relative to the basis of the holon and magnon states, 
the new basis will of course be incomplete. Special care 
will be taken to preserve the magnon states and their 
interaction with the holons. Only the first term on the 


right-hand side of Eq. (52) will be affected by the projec¬ 


tion, which comprises three steps. 


1. Logarithmic discretization of the conduction band 

The first step is controlled by an arbitrary dimension¬ 
less parameter A > 1 . Given A, we split the conduction 
band into two segments, one above and the other below 
the Fermi level, and divide each segment into an infinite 
logarithmic sequence of intervals 

l m ±=±D( A- m ~\A- m ) (m = 0,1,...), (57) 





















where the positive segment runs from e/~ = ep = 0 to 
6 k = D, and the negative one from = -D to = 0. 

For each interval , we define the normalized Fermi 
operator 

\rn/2 _ 

Qm± = —7= E ^ k ( m = 0, 1, - - ( 58 ) 

%/1 “ A keim± 

It follows that the operator T(0) is a linear combina¬ 
tion of the a m ±: 


* (0) = Vl-A-1 E A _m / 2 (a m+ + a m _) = \/2/ 0 , 

(59) 


m=0 


where we have defined the normalized Fermi operator / 0 , 
so that Eq. ([52]) can be written in the exact form 


H = E*I(i> - 2 uS z M T*)*k + j± (/oV/o^M + H. c. 

k 

+ 2J|| ^/o^/o^ — foifoij 
+ 2E ^/o^/ot + foifoi) • 


(60) 


Here t= R and L refer to the chirality pseudospins 
defined in Sec. |IV[ From now on we shall call these simply 
“holon spins”, whereas S m is the “magnon spin”. 

Next, we project the first term on the right-hand side 
of Eq. (60) on the basis of the a m± , which yields the 


approximate equality 


H = v(a ] 


-2 u(af m +a„ 










where 


mf mt 

+ J± (/of/oi^M + H - C ‘) 

+2-/|| (/of/ot - fUn) S'. 
+2V (/ 0 V/ ot + /oj,/oj,) 


£ - 1_ — _A _m 


M 


log A 


(61) 


(62) 


is the average momentum in the interval 

The projection on the incomplete basis intro¬ 
duces deviations in computed physical properties of 

7 f2 

0[exp(— -—— )]P 7 which are insignificant for A < 3. 


we define a new infinite basis {/ n } (n = 0,1,...), where 
the Fermi operators f n form an orthonormal sequence of 
linear combinations of the a m± , 


fn = E 
m =0 


^nm 


^ra+ T ( 1) 


(63) 


In the new basis, f Q is the operator defined by Eq. (59), 
so that 


aom = (64) 

and the remaining Fermi operators are tailored to the 
requirement that 


oo 


oo 


f m (a| n+ a m+ d\ n _CL rn _) — t n (fnf n 44 c.). 

m=0 n=0 

(65) 


Wilson has shown that Eq. (65) is satisfied witlfEO 


y(i - A _2n_i )(i - a- 2 "- 3 ) 


( 66 ) 


and determined the coefficients a nm . In particular, for 
n — 1, 


^lm — 


1 — A 


-3 


3m/2 


(67) 


for 
, to 

obtain the Lanczos transformed Hamiltonian 


We now substitute the right-hand side of Eq. 
the kinetic terms on the right-hand side of Eq. 


(65 

(61 


£=E*« 

n =0 


n'l'f n+lt "h ^n\,$ n+lt^ 


~ filfn+n) SZ M + H. C. 
+ '^-L C.) 

+ 2J|| (/ 0 Vo t — foifoij S Z M 
+ (/ot^ot + foifoi) ’ 


( 68 ) 


which leaves us in position to truncate the infinite basis. 


5. Truncation and scaling 


2. Lanczos transformation 


Numerical treatment of the Hamiltonian (61) calls for 


truncation of the infinite basis. To preserve the inter¬ 
action terms on the right-hand side, before truncation 


Inspection of Eq. (66) shows that the fraction on the 


right-hand side rapidly approaches unity as n grows, so 
that to an excellent approximation, 


t n = EqK n / 2 (n such that A n < 1). 


(69) 

























If we are interested in physical properties at the en¬ 
ergy scale E , it is safe to truncate the infinite series on 
the right-hand side of Eq. ( 68 ) at n = TV, where N is 


the smallest integer satisfying tw < cirE for a specified 
dimensionless infrared cutoff e IR <C 1. Given an energy 
E and an infrared cutoff e 7i? , we define the truncated, 
scaled Hamiltonian 


H n = 


1 


E>n -i 


' N—l 


n =0 


rt 


fn 


nf-Zn+lf 


+ flifn+li) 


H“ J_ 


Mfltfn+lt 

(/of/( 

2 / (+/, 


Ot^M 


ot 


flfn +U )S Z M + H- C. 

+ H. c.) 

/oVoi) 



(70) 


where V n denotes the right-hand side of Eq. (69), i. e., 
V n = f 0 A-"/ 2 . 


4- Iterative diagonalization 


The truncated form (70) is convenient for iterative 
diagonalization . 27 With N = 1, only the magnon spin 


and the operators / 0 and f 1 contributing to the right- 
hand side, the Hamiltonian Hn=i is equivalent to a 
32 x 32 matrix that is easily diagonalized, numerically. 
The next Hamiltonian in the iterative procedure, i 7 /v =2 is 
then projected on the basis | r)i, /t +2 f I r )i> /jv+ 24 . I r )i, 

/jv+ 2 t/iv+ 2 i I r )i ( r = 1,2,...,32), where | r)i is one 
of the eigenstates resulting from the diagonalization of 
Hn= 1 , and the resulting 128 x 128 matrix is numeri¬ 
cally diagonalized. This completes the first iterative cycle 
(AT = 2). 

More generally, given the eigenstates | r) n of i7/v, the 
Hamiltonian H^+i is projected on the basis resulting 

from the operators 1 , /t+i t , /jv+u> and /iv+it/t+U a P‘ 
plied to the | r) n and diagonalized. The number of eigen¬ 
states of Hn in the first few iterations is N r = 2 1+2Ar . 
To keep the rapidly growing N r under control, only the 
eigenstates of Hjsr corresponding to eigenvalues below 
a fixed ultraviolet cutoff E uv are computed at each it¬ 
eration, so that the computational cost grows linearly, 
instead of exponentially, with the number of iterations. 
The ultraviolet cutoff E uv , the discretization parameter 
A and the infrared cutoff e IR control the computational 
cost and the accuracy of the physical properties calcu¬ 
lated from the eigenstates and eigenvalues of Hjsf. 


5. Renormalization-group transformation 


The factor l/2£v-i on the right-hand side of Eq. (70) 
expresses the Hamiltonian Hn in units of its small¬ 


est matrix element, t^-i ~ P/v-i and defines a 
renormalization-group transformation T, which turns a 
scaled truncated Hamiltonian into another scaled Hamil¬ 
tonian, one that is truncated at a finer scale. More specif¬ 
ically, from Eq. (70) we have that 


H n+2 =A H 


N ' 


N+l 

£ 

n=N 


tn 


E>n+i - 


v (fUf*+ It + flifn+li) 


Mfl t fn+lt ~ flifn+u) S M + H- C. 


(71) 


which defines the transformation T[i7/v] = i7/v+2- 

At first sight, it may seem more natural to step N -+ 
N + 1. Nonetheless, the Hrt -+ i7/v+i transformation is 
unwieldy, because successive applications almost invari¬ 
ably generate two-point limit cycles.^ The Hjy i7/v+ 2 
transformation, by contrast, has simple fixed points. 


B. Fixed points 


As first shown by Wilson,^ the spin-flip term on the 
right-hand side of Eq. ( [52] ) is a marginally relevant oper¬ 
ator. To identify the fixed points of T[Hn] we will there¬ 
fore consider two forms of the Hamiltonian that remain 
invariant under scaling of J ±. For weak coupling between 
the magnon and the holons, for small N the scaled trun¬ 
cated Hamiltonian Hjy in Eq. (70) tends to be close to 


one of them. As N grows, our numerical analysis shows 
that Hn flows towards the second fixed point. We will 
therefore refer to the former (latter) as the weak-coupling 
(strong-coupling) fixed points. 


1. Weak-coupling fixed points (J± = 0) 


Without the spin-flip term, the Hamiltonian (70) reads 


N—l 


n =0 


Vn-\H% — t n r(/,t/ || + 1 * + fnlfn+ll) 


-Mfl t fn +It ~ flifn+u) S M + H. 

+ 2 J\\ (A/of — /of/oj.) ,C, 'm 

+ 2V (A/ot + /04./04.) ’ 


c. 

(72) 


a form that commutes with S^. 

The spectrum of H ^ can therefore be classified by the 
z component S M of the magnon spin. It can be divided, 
in other words, into an s z M = 1/2 sector and another with 

S M “ —1/2. 

For each magnon-spin component, the Hamiltonian in 
Eq. ( [72] ) splits into an f and a f terms, labeled by the 
holon spin. For s z M = +1/2, we have that 


K = H °+1 + H n + .u 


ME 


(73) 

















where 


N—l 


Vn-iH *+ =(v - u) E t n (fl t fn+lt + H. C.) 


n =0 


+ (2V + J||)/oVot 


(74) 


and 


N—l 


=(v + «) E tnifllfn+U + H - C ') 


n =0 


+ (2V«J|,)/U 




(75) 


With = —1/2, the Hamiltonian (73) likewise splits 
in two terms: 


rrO _ rrO— . rrO— 

— ^AT,t 


(76) 


where 


N—l 


V N - L H° na =(v + u) E ^(/If/n+lf + H. c.) 


n=0 


+ (2V - Jn)/ 0 Vt 


Ot-'Of 


(77) 


and 


N—l 


=(v - u) E tn(flf n+u + H. c.) 


n=0 


+ (2V + 




(78) 


Comparison of the right-hand sides of Eqs. (74) and ( |78| ) 
shows that and transform into each other un¬ 


der the inversion f n ^ CA f n ^ and therefore have the same 

spectrum. Likewise, and transform into each 

other under the same inversion and have identical spec¬ 
tra. 

In fact, the four Hamiltonians, and have 

the same general form 


H 


D,W 


N 


N—l 

Dj2tn(flfn + 1+K- C.) + WJqJq 
n =0 


(79) 


where D and W are constants. 

It has long been established that, as N grows so that 
N —>> N + 2, the Hamiltonian Hp w rapidly approaches 
a single-particle fixed point labeled by a phase shift S\ 47 ^ 
Although the fixed-point Hamiltonians depends on the 
parity of AT, the structures for even and for odd N are 
similar. If N is even, for example, the fixed-point Hamil¬ 
tonian is given by the expression 


m = D 


oo 

voglao + E ln™+9m+9m+ 

m= 1 



(80) 


where r]o is a number between — y/X and y/X that has to 
be determined numerically, and to a good approximation 
(within 0 ( A _2m ) relative deviation), 

Vm±= (m = 1,2,...) (81) 


with a phase shift S dependent on the ratio W/D : 


r A W 

tan o = —7r——. 

2D 


(82) 


Equations (80), (81), and (|82j) show that, as the even 


integer N grows to infinity, the scaled truncated Hamil¬ 
tonian Htf flows to a fixed point comprising two infinite 
sets of eigenvalues {s^ n ±} and {s^ n ±} (m = 0,1, 2,...) as¬ 
sociated with holon spins that are parallel or antiparallel 
to the magnon spin, respectively. 

The eigenvalues £^± form the infinite logarithmic se¬ 
quence 


4± = («-«) (83) 

where 

V±J\\/2 

tan5± = —7r--— (84) 

v^fu 


Since the phase shifts 5± can take any value from 
—7 t/2 to 7 t/2 , depending on V and Jj_, Eq. (80) defines 
a line of (weak-coupling) fixed points. Physically, the 
two sequences corresponds to the two bands discussed in 
Sec. |IV[ with velocities v d= u in the magnon reference 
frame. Each band is phase shifted. The phase shift is 5+ 
for holon spins that are parallel to the magnon spin, and 
S- for spins that are antiparallel. Since c) + ^ 5- each 
flip of the magnon spin triggers two Anderson orthogo¬ 
nality catastrophes, one for each component of the holon 
spin. We will come back to this issue in Sec. YE Before 
that, however, we have to examine the other set of fixed 
points. 


2. Strong-coupling fixed points (J± —>• oo) 


As the coupling J± on the right-hand side of Eq. (70) 


grows, the magnon and the f Q orbital lock into a sin¬ 
glet, with energy below the other configurations of 
the magnon f Q -orbital pair. As J± —>> oo, the other con¬ 
figurations cannot contribute to the physical properties of 
the Hamiltonian. The degrees of freedom associated with 
both the magnon and / 0 are lost at this strong-coupling 
fixed point. At iteration N the basis is therefore reduced 
to the operators f (n = 1,..., AT, fi =t, |). The scaled 
truncated Hamiltonian can again be split into two decou¬ 
pled components, i. e., we have that 


TT OO _ TTOO 

n N — n N, t 


H 




(85) 












where both components are of the general form (79). 
Specifically, we have that 


N—2 


rroo 

n N ,n 


Doo E tMUn+lv + H. c.) + Woo/oV, 


O/Lt? 


n=0 


( 86 ) 


where fi =t, i and and Woo are constants that must 
be extracted from the iterative diagonalization of Hjy. 


The series on the right-hand side of Eq. (86) runs from 
0 to N — 2 because we have shifted the summation index 
n 7i+l. The coefficients Dqq and are independent 
of the spin index /q because H°°’* must have the symme¬ 
try of Hjy, which remains invariant under ^-inversion. It 
follows that the single-particle spectra of H and 
are both constituted by energies D 00 fj rn ± {m = 0,1,...), 
where 770 , ± are numbers between 0 and a/A, and 


Vm± 


= A m +5 : Ff (m = 1,2,...). 


(87) 


Here again, S is a phase shift, given by the expression 

Woo 


tan S = —7r- 


Do 


( 88 ) 


Again we have a line of fixed points. In contrast with 
the weak-coupling fixed points, however, each strong¬ 
coupling fixed points comprises a single, spin-degenerate 
band, with the phase shift 5. This is consistent with the 
expectation that the magnon velocity, which controls the 
difference between the velocities v ± u in the magnon 
frame, flows to zero as Jy, J±_ flow to strong coupling, as 
discussed in Sec. hyh 


C. Deviations from the weak-coupling fixed point 

1. Deviations 


Before discussing the behavior of the Hamiltonian 
Hat,o in the vicinity of the weak-coupling fixed point, 
brief recapitulation of NRG linearization seems appropri¬ 
ate. Given a fixed point ft*, the renormalization-group 
flow of a Hamiltonian ft tv in the vicinity of ft*, to linear 
order in the deviation between ft n and ft* is described by 
the approximate form 

V N -iSh N = ^ w p O p , (89) 

p= 1,2,... 


where the coefficients w p depend on SJin, and the O p are 
the eigenoperators of the renormalization-group trans¬ 
formation T[ft/v], i. e., each O p is an operator with the 
symmetry of ft tv satisfying the equality 


T[Op] = XpOp , 


In practice, the eigenoperators are constructed from 
an infinite sequence of Fermi operators (f n (n = 0,1,...) 
defined by the equality 


E< 


9m+ T ( 1) 9m — 


(91) 


where the g m ± are the operators that diagonalize i7|, 
defined by Eq. (80), and the coefficients a nrn are those in 
Eq. (63); for particle-hole symmetric ft*, the (j) n coincide 
with the f n . 


Each operator O p is a linear combination of products 
of 0 n ’s and 0^’s respecting the symmetry of ft* . The re¬ 
sulting eigenvalue is given by a simple rule: 27 46 Uach op¬ 
erator (\> n or with even (odd) index n contributes a fac¬ 
tor A - ^ n+1 )/ 4 (A - ( 2n+1 )/ 4 ) to the operators on the right- 
hand side of Eq. ( [89] ) . The index n of the <fi n operators 
therefore defines a hierarchy; the most relevant Op s are 


the bilinear forms q</> 0 /D_i compatible with the sym¬ 
metry of ft*, which have eigenvalue A = 1, next come the 
bilinear forms P_i(0o0 1 + H. c.), with A = A -1 / 2 , then 
D_i(0o02 + H. c.) and D_i {(j>\(j)i +H. c.), with A = A -1 , 
and so successively. 


2. Weak-coupling fixed point 

Given the symmetry of the weak-coupling fixed points, 
which conserve charge and Sf ot and remains invariant 
under ^-inversion, its most relevant eigenoperators are 

or = —E(E! +4^, • Sm + H. c.), (92) 

®2 = 1 ^t^ 0 t ~ |)^M) (93) 


and 


O 


wc 

3 


1 

TA/v-i 




(94) 


The operator 


1 




H. c.), by contrast, 


Vn-i 

is not an eigenoperator, since it is the linear combina¬ 
tion of Oi° e G™ 0 , which have distinct eigenvalues: as 
Sec. VE will show, the eigenvalue Ai of Of c depends on 
the parameters of the Hamiltonian, while A 2 = 1. 

The deviations of the Hamiltonian ift/v from for 
small N are therefore approximately described by the 
corrections due to the perturbation SH^ C defined by the 
equality 


5H% C = w 1 or + w 2 Or + w 3 OZ c , (95) 


with eigenvalue A p , at the fixed point. 


( 90 ) 


with coefficients Wi (i = 1, 2, 3) that depend on the cou¬ 
pling constant J± and on the fixed-point phase shifts 5±. 


















3. Lanczos operators 


Equations (92)-(94) relate the Oj (j = 1,2,3) to the 
Fermi operators (p^~X n = 0,1,...). As Eq. ( |9l| ) shows, 
the latter are linear combinations of the Lan czos op er- 
ators f n . In particular, as explained in Sec. VC l| in 


the vicinity of particle-hole symmetric fixed points, with 
phase shifts ( s ) = 0orc) = 7r/2, we have that <fi n = f n 
(n = 0,1, 2 ,...). Since the weak-coupling fixed point is 
particle-hole asymmetric, we instead have that 


fm ~ 7 


(96) 


with coefficients 7 m?n dependent on the fixed-point phase 
shifts. For instance, 70,0 ~ cos£, while 770 ~ sin 5. 

It follows from Eq. ([96]) that 


fofo = 7o,o 0o0o + [7o,o7o,i (0i0o + H. c.) 

+ other irrelevant terms], (97) 

which shows that (1 /X>tv-i)/o/o * s marginal. 

Analogous considerations apply to operators of the 
form (l/P7v-i)/m/ m +i (jn = 0,1 ,...), for which we have 
that 

/m/m+l = 7m,07m+l,0 00^0 

A (7m,07m+l,l 00^1 "h 7m,l7m+l,0 01^0 
+ other irrelevant terms). (98) 

Under particle-hole symmetry, 7 mj0 vanishes unless 
m = 0. It follows that the first term on the right-hand 


side of Eq. ( [98] ) is zero and that (l/V N -i )/m/? 


m+l 


is ir¬ 


relevant in the vicinity of a particle-hole symmetric fixed 
point. If the fixed point is particle-hole asymmetric, by 
contrast, 7 m ,o7m+i,o 7^ 0 and (l/V N _ 1 )f} n f m+1 will be 
marginal. 


D. Weak-coupling fixed point with u — 0 


As N grows, iLjy approaches a fixed point of the 
form (80). Consider, now, the effect of the perturbation 


5Hn, u = — 


2 u 


N-l 


T*n- 


1 ™_n 


m =0 


+H. c.]. 


( 100 ) 


Since Hn ,o is particle-h ole as ymmetric, every term on 
the right-hand side of Eq. (100) is marginal, as shown by 


Eq. (98). Given that the coefficients 7 m; o decay rapidly 
with m, the contribution of (/q^/i^ — 
matrix elements of 5Hn, u is substantially larger than 
those of the terms with m > 0. To a good approximation, 
therefore, we can rewrite Eq. (100) in the form 


SH 


N,u = 


2 uto 

V N -i - 


(/oV: 


Ot-'lt 


foifli 


)S Z M + H. c. 


( 101 ) 


Substitution of Eq. (101) on the right-hand side of 


Eq. (70) yields the following approximate expression for 


the scaled truncated Hamiltonian 


N-l 


J M 


V N-lH N = vJ2 *n(/nt/n+.l t + /«7«+U + H ' C ‘) 

n=0 

- 2^o (/oV/it - /oVu + H ‘ C ) S '. 

+ ^(/oVo7m+ H - c -) 

+ 2 J|| (. 7 /ot - /04,/oi) ^ 

+ 21 7 (/oVot + / 0 V 04 ) • 


; M 


( 102 ) 


Section |VF| presents numerical data resulti ng fro m the 
iterative diagonalization of the Hamiltonian (102). Pre¬ 


liminary to that, however, it seems appropriate to qual¬ 
itatively discuss the renormalization-group flow on the 
basis of perturbative results. 


For small parameter u, an alternative view of the scaled 
truncated Hamiltonian Hn proves valuable. Recall that 
this condition is satisfied in the “slow magnon” regime 
u <C u, as discussed in Sec. Hy Aj Under this condition, 
we can regard both the terms proportional to u and J_\_ 
on the right-hand side of Eq. (70) as perturbations. The 


unperturbed Hamiltonian is then given by the expression 


H% = 


f N-l 


Vn- 


\ n =0 


- ] E [ V (fltfn+lt + flifn-m) + H - 


1 M 


2 ^ll (/oV/ot foifox) ^ 

' ^ (/ot^ot + fox fox) 1 * ( 99 ) 


E. Qualitative discussion of the 
renormalization-group flow 


To visualize the flow of the Hamiltonian (102) in 


renormalization-group space, it is convenient to split Hn 
into an unperturbed term and a perturbation defined by a 
linear combina tion o f the eigenoperator s Q™ c (n = 0,1, 2) 
defined in Sec. V C Inspection of Eq. ( 102 | ) identifies on 
the right-hand side the spin-flip term, that is, the term 
proportional to Jj_, in the absence of which the Hamil¬ 
tonian Hn can be easily diagonalized, as Sec. nm ex¬ 
plained. The spin-flip term can be regarded as (propor¬ 
tional to) a part of the o perat or 0f c . In order to split 
the right-hand side of Eq. (102) into an unperturbed part 


and a perturbation proportional to Of we can choose ei- 


































ther 


where 


Vn-iSHn = J\\ 


(/oVo^m + H. c.) 


+2(/o-|-/ot flifoOSk 


(103) 


or 


V n ~iSH n = J\\ 


(/oVoi S m + H. c.) 
2(/o-|-/ot — /o|/o|)5m 


(104) 


as the perturbation. 


With / 0 —>>00’ the Hamiltonian in Eq. (103) is evi¬ 


dently prop ortio nal to 0f c . To show that t 
nian in Eq. (104) is also proportional to O± 


to let S^ 


re Hamilto- 
one only has 

—S^, a gauge transformation that switches 
the sign of the first term on the right-hand side with¬ 
out affecting the second. We therefore have two choices. 
To decide between them, we have to examine the unper¬ 
turbed part of the Hamiltoni an re sulting fro m eac h al¬ 
ternative. Comparison of Eq. (102) with Eqs. (103) and 
( |104| ) shows that the coefficient of the term proportional 
to (/ 0 f t / 0t — /q|/o|)^m i n unperturbed Hamiltonian 
stemming from Eq. ( |103| ) is smaller than the coefficient in 
the unperturbed Hamiltonian associated with Eq. (104). 


It results that SHn is more relevant th an SH n and, there¬ 
fore, that Eq. (|104|) prevails over Eq. (1103). 


The unperturbed Hamiltonian is, then, 


N—l 


Vn-iH° n = V £ t n (fUf n+n + flf n+ u + H. C.) 

n=0 

- 2^t 0 (/ 0 f t / lt - / 0 \/ u + H. c }jS z M 

+ 2(J|| + J±) (/o f t /ot — /oVot) 

+ 2E ^/o^/ot + /o\/ot) ’ (105) 

and the perturbation, 

V N -iSH N = 2J± J2 fojou V/iv • S M > (106) 


so that Hj\f = Hn o + SHjy. 

In view of Eqs. (92), after the irrelevant contribu¬ 


tions on the right-hand side of Eq. (97) are disregarded, 


Eq. (106) can be written in the form 


Vn-iSHn = 2J±lloOT 


(107) 


Like the Hamiltonian ([72]), Hn o is easily diagonalized. 
Since S Z M is conserved, its eigenvectors can be classified 
by the z component of the magnon spin. Given the eigen¬ 


value the right-hand side of Eq. (|105|) is quadratic 


and hence numerically diagonalizable. Alternatively, we 


can split the right-hand side of Eq. (105) into an and a 
terms: 


77/vo = T7/vo,t + 77jvoq, 


(108) 


N—l 


HnO,/j — v tnifnnh 


n+l/x 


H. c. 


n —0 


+ 


2^/oVo m -^o(/oVi m +H. c.) 


+ (J\\+J±) 0 z flU (e z =±l), (109) 

and we have introduced the shorthand 


6 Z = (: 2s z M )(2fi) 1 


( 110 ) 


so that 0 Z = 1 (6 Z = —1) when the magnon and holon 
spins are parallel (antiparallel). 

For each 6 Z1 we can now let f n —)> f n _ 1 (n = 

1,2, ...,7V) and f Q —)> Co to identify 77 ato/x with a 
resonant-level Hamiltonian, i. e., a spinless U = 0 An¬ 
derson Hamiltonian in which Co represents the impurity 
level and the remaining Lanczos operators f n represent 
the holons. In this picture, the impurity has energy 
2 V + 0 Z {J\\ — J±) and its coupling to the conduction 
band is to(v — u0 z ). It follows that, for each product 
(2s z m ){2h) = +1 or (2s z M )(2fi) = — 1 , the conduction 
levels form a conduction band with uniform phase shifts 
given by the expression 


tan 8(6 z ) = — tvv 


V + 0 Z (J\\ + 7j_)/2 

(v - 0 z u ) 2 


(0 Z = ± 1 ). 

(Ill) 


The weak-coupling fixed point comprises distinct phase 
shifts. For instance, if the magnon spin is (J|), the t 
holons have the phase shif ts in Eq. (Ill) with 0 Z = 1 
(0 Z = —1). Equation (111) is the analogue of Eq. ([84]) . 


While the latter determines the phase shifts of the weak- 
coupling fixed point for the scaled truncated Hamilto¬ 
nian (70), the former determines the phase shifts for the 


approximate Hamiltonian (102). With u = J± = 0, the 


two expressions for the phase shifts coincide. For small 
u/v the first-order expansion of Eq. (84) reads 


tan £± = —7r 


V±Jn/2 


U) 


( 112 ) 


while the first-order expansion of Eq. (Ill) for J± = 0 
yields 

tan^) = - A + M||/2 (l + 2fl,-V (113) 
v V vJ 


Compared with Eq. (112), Eq. (113) shows that the 


truncation leading to Eq. (102) leads to renormalization 
of the anomalous velocity u, so that u 2 u. Otherwise, 
to first order in u/v, the conduction-band phase shifts are 
identical. We will next show that the phase shifts S(0 Z ) 
control the renormalization-group flow of Hn to insure 
that the numerical results in Sec. EE computed with the 
appr oxim ate form of the scaled truncated Hamiltonian in 
Eq. (102), yield quantitatively reliable conclusions con¬ 


cerning the magnon-holon interactions. 






































1. Perturbative treatment of the deviations from the 
weak-coupling fixed point 


The perturbation (106) couples states with different 
s z M and hence breaks the degeneracy between s z M =ff, 
g =1 states and s z M =4, /a states. Consider for ex¬ 
ample the following degenerate eigenstates of the unper¬ 
turbed Hamiltonian: 




(114) 


and 


\tA) = ffl+T I fU), 


(115) 


where and if indicate the z component of the magnon 
spin, and g\ +ll creates an holon at the first level above 
Fermi level with spin component g. 

To compute the first-order correction to the energies 
of the two states, we have to diagonalize the matrix 


V N -i 6HW = 


— 2J±7o ? o 


\ia\ot 

(tA\or 


i, a, ft i or 
ft) a, ft i or 


4 ft) 


(116) 


where 70,0 = cos£(—), because the holon states g\ + and 
the magnon spin are antiparallel both for 4,4) and |ft, if). 

The two diagonal elements are identical. Straightfor¬ 
ward computation shows that 


<4ftiori4ft)w = <t,ft 1 or it,ft) 


N 


—A - 
2 


N-l 

2 


with the coefficient 0^01 defined by Eq. (64) 


(117) 


The two off-diagonal matrix elements are also identical. 
Each can be factored into a matrix elements with holon 
spin components g =t and another with g —\ r . We for 
instance have that 


a, it-1 or 1 tA) N = <tu 1 or u,^> 


in 


— I Sh+^of An 
x (4b +1 <l>oi9i+i 14b — )am 

(118) 

where | Fl^Oz) (m =t4> @z = ±) denotes the ground 
state of the g component of t he co nduction band with 
phase shift S(6 Z ), given by Eq. S 


The two factors on the right-hand side of Eq. (118) 


being equal, we only have to consider the first one. Since 
the states | £2^, +) and | £4, ~) ^ ave distinct phase shifts, 
the matrix element expresses an Anderson orthogonality 
catastrophe and decays with N following the power lav® 


(£ 4 ^, | |£" 4 , n — <a(^)A 


Al-5/TT) 2 


(119) 


where 5 = 5(— )—5(+), and the N -independent coefficient 
a (5) satisfies a(5 = 0) = aoi; for 5 4 0 the coefficient 
can only be computed numerically. 

The phase-shift difference 5 de pends on the Hamilto¬ 
nian parameters. From Eq. (Ill), it follows that 5 > 0 
for 


J _L + J\ I > 


4 uv 


-V. 


( 120 ) 


Equations (117) and (119), determine the two eigen¬ 


values of the perturbative matrix (|116|): 
Se± 


J A. 


-4 cos 2 S(-) ah±a 2 (S)A 4 4 2 4 
Vi L 


( 121 ) 


where the eigenvalues £e + and £e_ correspond to the 
triplet and singlet combinations of 4, ft) and 4? ft)? re " 
spectively. 


If the parameters are such that inequality (120) is sat 


isfied, S will be positive. Equation (|121| then shows 
that the splitting between the triplet and the singlet will 
grow as a small positive power of the inverse energy scale 
1/Xhv-i- The operator Of vc is therefore relevant and will 
drive the scaled truncated Hamiltonian away from the 
weak-coupling fixed point. 

T he re normalization-group evolution expressed by 
Eq. (121) differs only quantitatively from the evolution 


for the standard Kondo model. We therefore expect the 
Hamiltonian to cross over from the weak-coupling to the 
strong-coupling fixed points under the influence of 0f c . 
Nonetheless, since perturbation theory is only reliable 
in the vicinity of the fixed point and hence inadequate 
to describe the crossover, numerical treatment becomes 
necessary. 


F. Numerical results 

The numerical diagonalization of the approximate ex¬ 
pression for the scaled truncated Hamiltonian Hjy in 
Eq. ([102]) follows the procedure outlined in Sec. |VA4l 
We choose the discretization parameter A = 3 and start 
out with N = 1, i. e., with the basis of 32 many-body 
states defined by the magnon-spin operator and the 
Lanczos operators /o and /1. The conservations of charge 
and ^-component of the total spin, and the invariance 
of the Hamiltonian under z — z inversion reduce the 
projection of Hi onto this basis to a block-diagonal form 
that yields to easy numerical diagonalization. The result¬ 
ing eigenvalues and eigenvectors seed the iterative cycle, 
which is stopped at N = 49, when Pjv-i becomes smaller 
than 10 _10 Z2. 

Ultraviolet truncation starts at iteration N = 4, when 
the highest scaled energies first exceed the ultraviolet cut¬ 
off E uv = 20 and are discarded. Our discussion being 
focused on renormalization-group flows, not on the com¬ 
putation of physical properties, the infrared cutoff e IR 





















needs not be defined. The iterative procedure is very ef¬ 
ficient: with A = 3 and E uv = 20, in a standard desktop 
computer a complete run takes times of the order of 100 
seconds. 


1. Flow in renormalization-group space 


To describe the trajectory of the scaled truncated 
Hamiltonian in renormalization-group space, this section 
presents numerical data for the A dependence of illus¬ 
trative scaled many-body energies resulting from the it¬ 
erative diagonalization of for various choices of the 
model parameters u , Jj_, Jy, and V. All parameters are 
expressed in units of v. Given that T[Hn] transforms Hn 
into Hn+ 2 , we plot the energies as a function of even it¬ 
eration number A. With A = even the spectrum of the 
strong-coupling fixed point is split into A/2 positive and 
A/2 negative single-particle energies, which makes devia¬ 
tions from particle-hole symmetry more visible than with 
A = odd. 


For clarity, for each run we show the A dependence of 
the lowest scaled energies Ejy in the (q = l,s z = 1/2) 
and (—1,1/2) sectors of the Fock space, where the charge 
q is measured from the charge of the half-filled band, 
the energies are measured from the ground state, and s z 
denotes the z component of the total spin. At the weak- 
coupling and strong-coupling fixed points, the two states 
correspond to simple holon configurations, as illustrated 
by Figs. [7] and [8] 

Figure [7] shows the two many-body states at the weak- 
coupling particle-hole symmetric fixed point. Since the 
number of single-particle states is N + 1, one of the levels 
lies at e = 0 and has ground-state occupation no = 1. 
A second holon (no holon) occupies the e = 0 level in 
panel (a) [panel (b)], which represents the lowest-energy 
state in the (1,1/2) [(—1,1/2)] sector. The two states are 
degenerate with the ground state. 


Figure [8] shows the two states at the corresponding 
strong-coupling fixed point. The magnon spin now being 
locked into a singlet with a holon, the remaining A single¬ 
particle levels are symmetrically distributed around the 
Fermi level. The minimum-energy state in the (1,1/2) 
[(—1,1/2)] sector contains a single particle (hole) at the 
lowest (highest) level above (below) the Fermi level. The 
two states are degenerate, with scaled energy f)i+ = ir¬ 
relative to the ground-state energy. The energies fji± are 
approximately described by Eq. (81), with S = Ett/2. 
For A = 3, in particular, fji± = 0.8. 


The following sections describe the A dependence of 
the lowest many-body energies in the (±, 1/2) sectors for 
various model parameters. For subsequent reference, we 
start out with an example of particle-hole symmetry. 
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Figure 7. (color online) Lowest-energy eigenstates in two 
sectors of the even-A" particle-hole symmetric weak-coupling 
fixed-point Hamiltonian with Jy = u = V = 0. Panels 
(a) and (b) show the lowest many-body eigenstate in the 
(q — l,s =s 1/2) and (—1,1/2) sectors, respectively. The 
horizontal lines represent the single-particle energies result¬ 
ing from the diagonalization of Hat, the blue spheres mark 
the filled levels, and the symbol indicates the z component 
of the magnon spin. Measured from the ground state, the 
many-body energies of both states are zero. 



Figure 8 . (color online) Lowest-energy eigenstates in two 
sectors of the even-A particle-hole symmetric strong-coupling 
fixed-point Hamiltonian with u — V = 0. Panel (a) shows the 
lowest many-body eigenstate at the (q = l,s = 1 / 2 ) sector, 
and panel (b) shows the lowest many-body eigenstate of the 
(—1,1/2) sector. The horizontal lines are the single-particle 
energies resulting from the diagonalization of Hn, and the 
blue spheres indicate the filled levels. Measured from the 
ground state, in which all negative levels are filled, both states 
have scaled energy 771 « 0.8, for discretization parameter A = 
3. 


2. Particle-hole symmetric Hamiltonians (u — V = 0/ 


With u = V = 0, Equation (52) describes a particle- 


hole symmetric, anisotropic Kondo Hamiltonian, the 
renormalization-group trajectory of which has long been 
known to merge onto that of an isotropic Hamiltonian. 49 
We therefore expect Hjy to flow from the vicinity of 
the particle-hole weak-coupling fixed point to the strong¬ 
coupling fixed point. Figure [9] shows an example, with 
J _= 0.2t> and Jn = O.Olv. The particle-like eigen- 

























































states in the (1,1/2) sector are degenerate with the hole¬ 
like eigenstates in the (—1,1/2) sector. The open cir¬ 
cles depicting the lowest-energy eigenvalue of Hjy in the 
(1,1/2) sector therefore coincide with the + signs de¬ 
picting the lowest-energy eigenvalue in the (—1,1/2) sec¬ 
tor as the Hamiltonian crosses over from near the weak- 
coupling fixed point to the strong-coupling fixed point. 
For comparison, the crossover for the isotropic model 
with J± = J\\ = 0.2 v is also shown. The inset shows that 
the crossover is universal, a horizontal shift AN = —9.1 
of the anisotropic curve being sufficient to make the two 
plots coincide. 



Figure 9. (color online) Scaled lowest energies En as functions 
of iteration number N for sectors (1,1/2) (open circles) and 
(—1,1/2) (crosses). The green circle s and crosses represent 
eigenvalues of the Hamiltonian ( |102| ) with u — V — 0, J± = 
0.2 v, and Jn t= 0.01 v. As N grows, the renormalization- 
group transformation T[Hn] drives the Hamiltonian from the 
vicinity of the weak-coupling fixed point, where E% — 0, to 
the strong-coupling fixed point, where E% = 0.8. For compar¬ 
ison, the evolution of the eigenvalues for a symmetric Hamil¬ 
tonian, with J|| = J x — 0.2 v, is also shown. The inset shows 
that, horizontally shifted by AN = 9.1, the isotropic curve 
(red solid line) coincides with the anisotropic curve (green 
crosses and open circles). 


3. Particle-hole asymmetry 


Figure [lO] shows the N dependence of the minimum en¬ 
ergies in the (1,1/2) and (—1,1/2) sectors for J± = 0.2 v, 
J\\ = 0.01 v, and V = 0.10 with u = 0. T he te rm pro¬ 
portional to V on the right-hand side of Eq. ( |102| ) breaks 
the degeneracy between the (1,1/2) and (—1,1/2) ener¬ 
gies and shifts the phases of the fixed-point single-particle 


levels, as discussed in Section |VB| The phase shifts at 
the weak-coupling fixed point affect the renormalization- 
group flow and delay the crossover to the strong-coupling 
fixed point. Universality is nonetheless preserved, as indi¬ 
cated by the congruence in the inset, which compares the 
isotropic curve in Fig. [9] with the two particle-hole asym¬ 
metric curves, both shifted horizontally by AN = —11.7 
and displaced vertically to insure agreement at large N . 



Figure 10. (color online) Scaled lowest energies in sectors 
(±1,1/2) as a function of iteration number for the indicated 
Hamiltonian parameters. As in Fig. [9] the corresponding en¬ 
ergies for the isotropic, particle-hole symmetric Hamiltonian 
are shown for comparison. The symbol convention follows 
Fig. IH The scattering potential V breaks particle-hole sym¬ 
metry and shifts the (1,1/2) [(—1,1/2)] energies up (down). 
The crossover between the weak-coupling and the strong¬ 
coupling fixed points nonetheless follows the universal be¬ 
havior in Fig. [9] as evidenced by the inset, which show the 
particle-hole symmetric, isotropic curve (red, solid line), hori¬ 
zontally displaced by AN = 11.7, and the asymmetric curves 
(green open circles and crosses) shifted vertically to make 
the large N energies match the particle-hole symmetric fixed- 
point energies. 


Figure [TT] shows results for u 7^ 0. The minimum en¬ 
ergies in the (±1,1/2) sectors for J_\_ = 0.2 v, J\\ = 0.2 v, 
and V = 0.10 v, with u = 0.2 v are plotted as functions 
of the iteration number N. Comparison with the flow in 
Fig. [lO] show that the spin-dependent contribution to the 
bandwibandwidthdth, i. e., t he te rm proportional to u 
on the right-hand side of Eq. (102) perturbs two features 
of the renormalization-group flow: the weak-coupling 
(strong-coupling) phase shifts S wc (S sc ) is closer to its 
isotropic value S wc = 0 ( S sc = 7t/ 2), and the crossover 
from the weak-coupling to the strong-coupling fixed point 
is shifted to the right. The delayed crossover indicates 
that the operator 0f c , defined in Se c. |V C , is less rele¬ 
vant. This conclusion agrees with Eq. (111), which shows 
that u reduces S(-) and increases <S(+) and therefore di¬ 
minishes the difference 6 , whic h co ntrols the exponent 
on the right-hand side of Eq. dm} . The ^-dependent 
ground-state phase shifts and the universal crossover con¬ 
firm that the term proportional to u in the Hamiltonian is 
a marginal operator, which displaces the high-energy and 
the ground-state Hamiltonians along the lines of weak- 
and strong-coupling fixed points, respectively. 

Analogous results are obtained for other parametrical 
choices. In particular, as u grows at fixed J_ j_, Jy, and 
V, the crossover is displaced to progressively higher N, 
while the strong-coupling phase shifts S sc are slightly dis¬ 
placed to higher or lower values, chiefly depending on V. 
Figure 12 shows S sc as a function of the ratio u/v for 
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Figure 11. (color online) Scaled minimum energies in the 
sectors (±1,1/2) as functions of the iteration number for 
the indicated Hamiltonian parameters. We follow the sym¬ 
bol convention in Fig. [9] and for comparison reproduce the 
isotropic curve in Fig. [9] As in Fig. [To] the energies in the 
(1,1/2) [(—1,1/2)] sector are shifted up (down) relative to the 
particle-hole symmetric energies. To show that universality is 
preserved, the inset shows the particle-hole symmetric curve 
(red, solid line) displaced horizontally by AN =12, superim¬ 
posed upon the particle-hole asymmetric curves (green, open 
circles and crosses) displaced vertically to insure agreement 
at the large N limit. 


fixed J_\_ = 0.2?; and illustrative choices of Jy and V. For 
u = 0, we expec t the weak-coupling phase shifts S(6 Z ) to 
be given by Eq. (121), and from the Friedel sum rule, ex¬ 


pect the strong-coupling phase shift to differ from S(0 Z 
by 7 t/2 . We therefore expect S sc to obey the relation 


V 

Cot &sq — 7T , 

V 


( 122 ) 


which is approximately satisfied for all data in Fig. [12] 
As u grows, 5 sc diminishes for small V and grows as 
V exceeds 0.05?;. In all cases, however, S sc depends 
weakly on u and is nearly independent of either J± or Jy. 
Roughly speaking, therefore, the strong-c oupli ng phase 
shifts are approximately described by Eq. (122) over the 
parametrical space of the Hamiltonian. Although the 
approximation distinguishing Eq. (102) from the scaled 


truncated Hamiltonian (70) restrict our analysis to small 


u, these results ratify the conclusion that the anoma¬ 
lous velocity is a marginal operator whose effect upon the 
physical properties of the system under study is limited 
to reducing the Kondo temperature and (slightly) dis¬ 
placing the low-temperature Hamiltonian along the line 
of strong-coupling fixed points. 


VI. CONCLUSIONS 

In this work, we considered the physics of doping 
a one-dimensional ferromagnetic insulator. Assuming 
spin-charge separation in the strongly interacting regime, 


Figure 12. (color online) Strong-coupling fixed-point phase 
shifts as a function of the anomalous velocity u at fixed J± = 
0.2?;, Jy = 0.1 v, and V = 0.1 v. The phase shifts are obtained 
from Eq. (87), the single-particle energies rjm± being identified 
as the limit to which the eigenvalue of the many-body state g ^ 
(m = 2) —sector (1,1/2) tends as N grows past the crossover. 
The isotropic limit of the phase shift, not shown, is S/tt = 0.5. 


we modelled the charge excitations by spinless fermions 
(holons) and the spin excitations by magnons. 

The interaction between magnons and gapless holons 
in the limit of small doping gives rise to infrared singulari¬ 
ties in perturbation theory. The mechanism behind these 
singularities is clarified by mapping to an effective field 
theory equivalent to an anisotropic Kondo model. In our 
case, the holons play the role of the conduction electrons 
and a magnon with momentum commensurate with the 
holon Fermi momentum pp plays the role of a mobile im¬ 
purity. The internal degree of freedom of the impurity is 
a chirality pseudospin, which indicates whether the mo¬ 
mentum is closest to pf or — pp. A pseudospin flip cor¬ 
responds to a 2 pp scattering process. Another difference 
from the conventional Kondo model is that the impurity 
(magnon) is mobile rather than localized. By applying 
a Galilean transformation to the frame moving with the 
magnon velocity ??, we find that the effective Hamilto¬ 
nian contains two holon bands with velocities v ± u or 
v — u depending on whether the holons move in the same 
or opposite direction as the magnon. 

We investigated the effects of the magnon-holon inter¬ 
action in the nonperturbative, low-energy regime using 
the numerical renormalization group. We generalized the 
numerical method to treat the effects of a new marginal 
perturbation proportional to the magnon velocity u. The 
renormalization group flow takes the Hamiltonian from 
a line of weak-coupling fixed points characterized by two 
distinct phase shifts (related to scattering between the 
moving magnon and the two holon bands) to another 
line of strong-coupling fixed points with a single phase 
shift 5 SC . The existence of only one phase shift can be 
interpreted in terms of the vanishing group velocity of 
the single-magnon excitation when the magnon momen¬ 
tum approaches pp. By analogy with the original Kondo 






















effect, we expect that at the strong-coupling fixed point 
the magnon pairs up with a holon to form a singlet of the 
chirality pseudospin. On the other hand, we should stress 
that, since the total magnetization is conserved, the real 
spin carried by the magnon is not screened. The resid¬ 
ual potential scattering between the remaining holons 
and the magnon-holon pseudospin singlet accounts for 
the nonuniversal phase shift at the strong-coupling fixed 
point. Finally, we showed that the initial value of u at 
the weak-coupling fixed point has a similar effect to the 
potential scattering term, namely it drives the phase shift 
away from the particle-hole symmetric value S sc = tv/2. 

The results of this work suggest that it would be inter¬ 
esting to study dynamical response functions probing the 
spin excitation spectrum of one-dimensional itinerant fer- 
romagnets when the magnon momentum approaches the 
holon Fermi momentum. In particular, the single phase 


shift 5 SC should have consequences for the power-law sin¬ 
gularity at the single-magnon threshold in the trans¬ 
verse spin structure factor, which is equivalent to the 
magnon spectral function. 23 Another interesting question 
is whether the Kondo effect discussed here can drive a 
nontrivial instability of the metallic ferromagnetic phase 
of a lattice model such as Tasaki’s model in Eq. ([!]), since 
it is expected to lower the energy of the single-magnon 
excitation. The behavior of the magnon spectral func¬ 
tion in lattice models, such as Tasaki’s model in Eq. 0 . 
could be investigated numerically using time-dependent 
density matrix renormalization group methods: 50 51 
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